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Abstract. - We study antiferromagnetism and single-particle properties in the two-dimensional 
half-filled Hubbard model at low temperature. Collective spin fluctuations are governed by a 
non-linear sigma model that we derive from the Hubbard model for any value of the Coulomb 
repulsion. As the Coulomb repulsion increases, the ground state progressively evolves from a 
Slater to a Mott-Heisenberg antiferromagnet. At flnite temperature, we flnd a metal-insulator 
transition between a pseudogap phase at weak coupling and a Mott-Hubbard insulator at strong 
coupling. 



Introduction. The Hubbard model and its generalizations play a key role in the descrip- 
tion of strongly correlated fermion systems such as high-Tc superconductors, heavy fermions 
systems, or organic conductors p]. Despite its simplicity (the model is defined by two pa- 
rameters, the inter-site hopping amplitude t and the local Coulomb interaction U, and the 
symmetry of the lattice) , exact solutions or well-controlled approximations exist only in a few 
special cases like in one-dimension (ID) or in the limit of infinite dimension P]. 

It is now well established that the ground state of the 2D half-filled Hubbard model on a 
square lattice has antiferromagnetic (AF) long-range order. The origin of antiferromagnetism 
is believed to depend on the strength of the Coulomb repulsion. At weak coupling {U ^ i), 
a Fermi surface instability gives rise to an insulating spin-density-wave ground state as first 
suggested by Slater [2]. In 2D, thermal (classical) fiuctuations preclude a finite-temperature 
transition (Mermin- Wagner theorem) and the phase transition occurs at Tn = 0. Since the 
fermion spectral function is gapped at T = 0, one expects, from a continuity argument, that 
it will exhibit a pseudogap at finite temperature as a result of strong AF fiuctuations [3] . At 
strong coupling {U ^ t), the system becomes a Mott-Hubbard insulator below a temperature 
of the order of U. The resulting local moments develop AF short-range order at a much 
smaller temperature J — it'^ /U), and AF long-range order sets in at Tn = 0. 

Although the Hartree-Fock (HF) theory gives a reasonable description of the AF ground- 
state, it fails in 2D since it predicts AF long-range order. Several alternative approaches, 
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which do satisfy the Mermin- Wagner theorem, have been proposed: Moriya's self-consistent- 
renormaUzed theory the fluctuation exchange approximation (FLEX) 0, and the two- 
particle self-consistent approach However, these approaches are restricted to the weak 
to intermediate coupling regime {U ^ At). The strong-coupHng regime is usually understood 
from the Heisenberg model for which various methods are available JJ. 

In this Letter, we describe a theoretical approach which provides a unified view of an- 
tiferromagnetism and single-particle properties in the 2D half-filled Hubbard model at low 
temperature (including T = 0) and for any value of the Coulomb repulsion. It is based on 
a non-linear sigma model (NLctM) description of spin fiuctuations. Since it takes into ac- 
count only directional fiuctuations of the AF order parameter, it is valid below a crossover 
temperature Tx which marks the onset of AF short-range order. 

Besides its validity both at weak and strong coupling, our approach differs from previous 
weak-coupling theories 0013] by the fact that it is a low-temperature expansion (0 < T ^ 
Tx). In particular, the fermion spectral function is obtained from a spin-rotation-invariant 
perturbative expansion around the (gapped) HF ordered state. This should be contrasted 
with perturbative treatments applied to (gapless) free fermions interacting with soft collective 
spin fiuctuations where no small expansion parameter is available 0. In Ref. jHI, one of the 
present authors reported a calculation of the spectral function in the weak coupling limit of 
the Hubbard model using a NLcrM description of spin fiuctuations. However, the limitations 
encountered by previous weak-coupling theories could not be fully overcome. 

Non-linear sigma model. The Hubbard model is defined by the Hamiltonian 

H = -t {cl,Cr'a+h.C.) + Uj2nr^nrl- (1) 

{r,r'),cr r 

cj^ (crcr) creates (annihilates) a fermion of spin a at the lattice site r. riro- = cJ^Crtj and (r, r') 
denotes nearest-neighbor sites. We take the lattice spacing equal to unity and ?i = /cb = 1. 

We express charge and spin fiuctuations in terms of auxiliary fields. For this purpose, we 
write the interaction term in iQJ as nr^nrx = [(cjcr)^ — (cJcr-r2rCr)^]/4 where fir is an arbitrary 
unit vector JHI- Cr = (crf , Crx)"^ and cr = (cti, (J2, 0-3) denotes the Pauli matrices. Spin-rotation 
invariance is made explicit by averaging the partition function over all directions of fir- In 
a path-integral formalism, fir becomes a time-dependent field. Decoupling the interaction 
term by means of two real auxiliary fields, Ac and As, the partition function is then given by 
Z = J 2?[ct,c]/ ^[Ac, As,rJr]e-^ with the action (/3 = 1/T) 



S — Si 



'0+/ dr^ 



A2 



U 



(2) 



where Cr, cl are Grassmann variables. Sq is the action in the absence of interaction. Eq. ^ 
defines an "amplitude-direction" representation, where the AF order parameter field is given 
by Asr^^r- In the following, charge fiuctuations (Ac) are considered at the saddle-point (i.e. 
HF) level: -lAcr = {U/2) (cjcr) = U/2. Below the HF transition temperature T^^ , the 
amplitude of the order parameter takes a well-defined value so that we can consider Asr 
within a saddle-point approximation, i.e. A^r = Aq where the fiuctuations of Ao are ignored. 
When T < T^f, Aq ~ te-^'^v^ for J7 < < and tends to U/2 for f7 > i. Below the 
crossover temperature Tx (to be defined more precisely later) which marks the onset of AF 
short-range order, the fir field can be parametrized by fir = (— l)'"nr(l — L^)^/^ + Lr, where 
the Neel field rir (|nr| = 1) is assumed to be slowly varying The small canting vector 

Lr, orthogonal to rir, takes account of local ferromagnetic fiuctuations. It turns out to be 
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Fig. 1 - Spin- wave velocity c, bare spin stiffness Ps, and fraction no of condensed bosons at T = 0. 



convenient to introduce a pseudo-fermion 0, 
Neel field, (pr — {(pri , 'Pri)'^ is defined by Cr 
SU(2)/U(1) matrix satisfying Rr<J-iR\ = cr ■ 



a whose spin is quantized along the (fiuctuating) 
= i?r0r where i?r is a site- and time-dependent 
riv. The action can then be written as 



r dTy24{dr-Aor-2t V COsHa^- V)-Ao[(- 1)^3^1^ +lr-'T] I ^ 



,(3) 



where we have introduced the SU(2) gauge field Aqt = —RldrRr, = iR^d^Rr {pi — x,y), 
and the rotated canting field Ir — TZ~^Lr. Here TZr is the S0(3) element associated to Rr 
which maps z onto rir. In Eq. (pj, both 1 and A^j, are small, since the gauge field is of order 
dfj^n. The effective action of the Neel field is obtained by expanding ijSJ to second order in 
these variables and integrating out the fermions and the canting field Ir. Skipping technical 
details, we obtain 



•Snlo-m [n] 



2 



dT 



(Vrir 



(4) 



where we have taken the continuum limit in real space. The bare spin stiffness and the spin- 
wave velocity c are given by = Cc/S and = (ec/2)(x2^ — U/2), where Ec is the absolute 
value of the (negative) kinetic energy per site and x± the transverse spin susceptibility in 
the HF ground state (Fig. In the weak-coupling limit, AF short-range order cannot be 
defined at length scales smaller than ~ t/^o which corresponds to the size of bound 
particle-hole pairs in the HF ground state. Thus Eq. Q should be supplement with a cutoff 
A ^ min(l,^(7^) in momentum space. The NLcrM I^J was first obtained by Schulz 22]. In the 
Hmit U ^ t, it reproduces the result obtained from the Heisenberg model with an exchange 
coupling J = At'^/U 0. 

Magnetic phase diagram. We solve the NLcM within the CP^ representation, where the 
Neel field is expressed in terms of two Schwinger bosons: — z^cTZr {zr = {z^^ , z^i)'^) 
with z^Zr = 1. When the CP^ representation is generaHzed to the CP^^^ representation by 
introducing N bosons Zra (o" = 1, ...,iV), the NLtrM can be solved exactly by a saddle-point 
approximation in the N ~* oo limit [Jj. At zero temperature, there is a quantum critical point 
at g = (7c = 47r/A between a phase with AF long-range order and a (quantum) disordered 
phase, g = c/p'^ is the coupling constant of the NLcrM. In the ordered phase {g < gc), 
a fraction uq — 1 — g/gc (0 < no < 1) of the bosons condenses in the mode q = 0. no 
determines the mean- value of the Neel field: |(nr)| = no. We find that there is AF long- 
range order for any value of the Coulomb repulsion in the ground state of the 2D half-filled 
Hubbard model. At weak coupHng, 1 — no = (?/ffc ~ e"^'^^^*^ is exponentially small. By an 
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Fig. 2 - Phase diagram of the 2D half-filled Hubbard model. All lines, except Tn = (thick solid 
line), are crossover lines. The NLcrM description is valid below Tx. FL - Fermi liquid phase, PG - 
pseudogap phase. The vertical dotted line indicates the finite-temperature metal-insulator transition 
obtained from the criterion p{ijj = 0) = 0. 



appropriate choice of the cutoff A, we reproduce the result uq ~ 0.6 for C/ » i as obtained 
from the 2D Heisenberg model on a square lattice pHI (Fig. GJ. At finite temperature, the 
AF long-range order is suppressed {uq — 0). However, the AF correlation length remains 
exponentially large: ^ ~ {c/T)e^'"''=^'^ , where ps = p^{l — g/gc) is the zero temperature spin 
stiffness. This regime, which is dominated by classical (thermal) fluctuations (since c/^ -C T), 
is known as the renormalized classical regime. The NLcrM description is valid for T < Tx 
when amplitude fluctuations of the AF order parameter are frozen and the assumption of 
AF short-range order holds (i.e. ^ ^ A~^). At weak coupling Tx ~ , while Tx ~ J at 
strong coupling. The phase diagram is shown in Fig. [2 Above T^^ , spin fluctuations are 
not important and we expect a Fermi liquid (FL) behavior. Between T^^ and Tx (a regime 
which exists only in the strong-coupling limit), local moments form (^o ~ 1) but with no AF 
short-range order (Curie spins: ^ '--^ 1). Below Tx, the system enters a renormalized classical 
regime of spin fluctuations where the AF correlation length becomes exponentially large. AF 
long-range order sets in at Tn = 0. Although there is a smooth evolution of the magnetic 
properties as a function of U, the physics is quite different for U -^t and U ^ t. This is shown 
below by discussing the fermion spectral function. The main conclusions are summarized in 
Fig. El 

Fermion spectral function. Now we consider the effect of long- wavelength spin fluctuations 
on the fermion spectral function. The fermion Green's function Q{v, r; r', r') — —{cr{T)c\, (r')), 
written here as a 2 x 2 matrix in spin space, is computed using Cr — Rr(t>r with (-Rr)tT = 
(i?r)J| — and (i?r)iT = ~ -^""i- Integrating flrst the pseudo-fermions and the 

canting fleld Lr, we can write the Green's function as 

g{l, 2) - I y PWe-^--'Hi?ig(l, 2\z)rI (5) 

where we use the shorthand notation i = (r^jT^). Q{1, 2\z) is the pseudo-fermion propagator 
calculated for a given conflguration of the bosonic fleld z. The action (pj, when expanded 
in powers of Afj_ and Lr, can be written as S'hf[</>] + S"[2;,0, L]. It describes HF pseudo- 
fermions interacting with spin fluctuations via the action S' . Since the HF pseudo-fermions 
are gapped, we expect a perturbative expansion in S' to be well-behaved. To leading order 
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Fig. 3 - Left: spectral function A{k,w) in the weak-coupling limit U = t for T — (Slater an- 
tiferromagnet) and T = Ao/5 (pseudogap phase). Right: spectral function in the strong-coupling 
regime U = 12t for T = (Mott-Heisenberg antiferromagnet) and T = J/5 (Mott-Hubbard insula- 
tor), k = (7r/2,7r/2). For T = 0, the vertical lines represent Dirac peaks of weight no/2. Note the 
difference in the energy scale, which is fixed by Ao, between the two figures. 



5(1, 2|z) g"P(l,2) with the HF Green's function. From I©, we then obtain 

5.(k,k',c.) = / t/HF(k-q,k-q,a;-c.,)P(q,c..) + noa?F(k,k',c.), (6) 

(k,k,a;) ^ -5,,,,-^ -,4,,,+^-^, p(q,^,)^„__, (7) 

where = /^^ ^ J^^ tt = (tt, tt), and a; (w^) denotes a fermionic (bosonic) Matsubara 
frequency. 2?(q, w^) is the Schwinger boson propagator (for q ^ 0) obtained from the saddle- 
point solution of the NLcrM. Here = c{a^ + ^-2/4)1/2 ^ (^2 _^ Ag)!/^, with 
Ck = —2t{coskx + cos ky) the energy of the free fermions. At finite temperature, uq vanishes 
so that the fermion Green's function is spin-rotation and translation invariant. From we 
obtain the spectral function .4(k, oj) = — 7r~^ImfJo^(k, \i,iuj u; + jO+): 

y4,(k,cj) y4,inc(k,Lj) + no^HF(k,t^), 

Anc(k,Cj) = y ^|[71_B(cJq) + ?lF(-£;k-q)] [ul_^5{uj - U;c^ - E]^^q) + vl_^5{uj + UJq + El^^ 

where npiijj) and nsito) are the usual Fermi and Bose occupation numbers and Ahf the HF 

spectral function: ^HF(k, w) = u^(5(a' — E^) + v^6{lu + E^), where u^,v^ = ^ ^1 ± • 

The normalization of the spectral function, J dLLjA{k,uj) = 1, follows from the saddle-point 
equation of the NLcrM. 

The spectral function is shown in Fig. El At weak coupling and zero temperature, our 
theory describes a Slater antiferromagnet. The AF gap 2Ao ~ te~'^^^*^^ is exponentially 
small. There are well-defined Bogoliubov quasiparticles (QP's) with energy ±i5k, as in the 
HF theory, but their spectral weight is reduced by a factor no due to quantum spin fiuctua- 
tions. The remaining weight (1 — no ~ g-'^^^/t/u ^ 2) is carried by an incoherent excitation 
background at higher energy {\uj\ > i?k)- As U increases, the Slater antiferromagnet progres- 
sively evolves into a Mott-Heisenberg antiferromagnet with a large AF gap and a significant 
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Fig. 4 - Density of states p(ui) at low energy. As the Coulomb repulsion U increases through a critical 
value Uc 2i 4.25t the pseudogap becomes a Mott-Hubbard gap [see the vertical dotted line in Fig. [5]. 

fraction of spectral weight transferred from the Bogohubov QP's to the incoherent excitation 
background. At strong coupling, the AF gap 2Ao ~ U and ~ 0.6. 

At finite temperature, A{k,Lu) exhibits two broadened peaks at the HF QP energy ±i?k- 
In the vicinity of the peak at i?k, the spectral function is well approximated by Ap^^^.{'k, uj) = 
w^^nsdw — i?k|)- The peak has a width of the order of the temperature and therefore corre- 
sponds to incoherent excitations. We find that the spectral weight of the peak, J dojAp^g^^Ck, uj), 
is temperature independent and equal to u'^m, which is nothing else but the Bogoliubov QP 
weight in the ground state. We conclude that the peak at i^k is an incoherent precursor of the 
zero-temperature Bogoliubov QP peak. As the temperature decreases, it retains its spectral 
weight but becomes sharper and sharper, and eventually becomes a Dirac peak at T = 0. As 
expected, the weak coupling pseudogap continuously evolves into the AF gap when T ^ 0. 

Eq. lO shows that the contribution to A{k, uj) at low energy involves the Bose occupation 
number ^^(wq). This indicates that the low-energy fermion states (0 < < i?k) are due to 
thermal bosons, i.e. thermally excited spin fluctuations. A fermion added to the system with 
momentum k and energy |a;| < Sk can propagate by absorbing a thermal boson of energy 
ojq and emitting a pseudo-fermion with energy i?k-q = + cjq. The lowest fermion energies 
are obtained by solving to = i?k-q — ojq (or lo = — i?k-i-q + Wq)- In the weak coupling limit 
(Fig. El, maxq(cjq) — cA ^ 2Ao and i?k-q ^ E\^- Thus there is spectral weight at zero 
energy: the spectral function and the density of states exhibit a pseudogap. Nevertheless, 
the density of states p{uj) — J^. -4(k, u) remains exponentially small at low energy: p{uj) ~ 
e"'^"/"^ cosh((jj/r), \uj\ <C Aq. This result differs from pseudogap theories based on Gaussian 
spin fluctuations which flnd a much weaker suppression of the density of states . It bears 
some similarities with the result obtained by Bartosch and Kopietz for fermions coupled to 
classical phase fluctuations in incommensurate Peierls chains In the low temperature 

regime dominated by directional fluctuations of the order parameter, the suppression of the 
density of states at low energy is indeed expected to be exponential. In the strong coupling 
limit, since ^ U/2 and cA ^ J <C J7/2, there is a gap (of order U/2) in the spectral function 
and the density of states. Thermally excited spin fluctuations reduce the zero-temperature 
AF gap U/2 by a small amount of the order of J. The system is a Mott-Hubbard insulator 



We therefore conclude that our approach predicts a flnite-temperature metal-insulator 
transition between a pseudogap phase and a Mott-Hubbard insulator as the strength of the 
Coulomb interaction increases (Fig. [JJ: at a critical value Uc, the density of states at zero 
energy p{u; = 0) vanishes and the pseudogap becomes a Mott-Hubbard gap (Fig. Uc is 
obtained by equating the minimum energy Aq of a HF fermion to the maximum energy of a 



(Fig. El. 
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Schwinger boson ^/m"^ + c^A^. For T ^ the result is Uc — 4.25i. However, being a low- 
energy theory, the NLcrM does not allow us to describe accurately the high-energy Schwinger 
bosons (with |q| ~ A) and in turn the low-energy fermion excitations. In particular, the 
critical value of U calculated above and the precise form of the density of states near zero 
energy, plotted in Fig. ^ depend on the cut-off procedure used in the NLctM. Note also that 
we do not know at which temperature and how the metal-insulator transition ends. 

Conclusion. We have presented a low-temperature approach to the 2D half-filled Hubbard 
model which allows us to study both collective spin fiuctuations and single-particle properties 
for any value of the Coulomb repulsion U. At zero temperature, it describes the evolution 
from a Slater to a Mott-Heisenberg antiferromagnet. At finite temperature, it predicts a 
metal-insulator transition between a pseudogap phase at weak coupling and a Mott-Hubbard 
insulator at strong coupling. Since the charge auxiliary field Ac is considered at the HF level, 
some aspects of the Mott-Hubbard localization are not taken into account. In particular, at 
intermediate coupling U ~ 8t, we expect both Bogoliubov QP bands (or precursors thereof at 
finite temperature) and Mott-Hubbard bands in the spectral function The Mott-Hubbard 
bands have a purely local origin, independent of the Fermi surface geometry, and should show 
up at a; ~ ±U /2 (with U/2 > Aq) in the spectral function. Nevertheless, we believe that our 
theory captures the main features of the physics of the 2D half-filled Hubbard model. 

On the basis of a numerical calculation (dynamical cluster approximation) , Moukouri and 
Jarrell have called into question the existence of a Slater mechanism in the 2D Hubbard 
model 51]. Using the criterion p{uj = 0) < 10~^/(2i) to identify the Mott insulating phase, 
they concluded that the system is always insulating at low (but finite) temperature even in 
the weak-coupling limit. From the same criterion (p(w = 0) < 10-2/(2^), 

we obtain a similar 

line in the (C/, T) plane as Moukouri and Jarrell. This shows that the numerical results of 
Ref. are not in contradiction with the existence of a Slater scenario at weak coupling, but 
refiect the exponential suppression of the density of states due to the presence of a pseudogap. 

* * * 

We thank A. -M.S. Tremblay for useful correspondence and S. Pairault for a critical reading 
of the manuscript. 

REFERENCES 

[1] For a general introduction to the Hubbard model, see F. Gebhard, The Mott metal-insulator 

transition (Springer- Verlag) 1997. 
[2] See, for instance, M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. Phys., 70 (1998) 1039. 
[3] J.C. Slater, Phys. Rev., 82 (1951) 538. 

[4] Y.M. ViLK and A. -M.S. Tremblay, J. Phys. I (France), 7 (1997) 1309. 

[5] T. MoRiYA AND K. Ueda, Rev. Mod. Phys., 49 (2000) 555. 

[6] N.E. Bickers and D.J. Scalapino, Ann. Phys. (N.Y.), 193 (1989) 206 

[7] A. Auerbach, Interacting Electrons and Quantum Magnetism (Springer- Verlag, New York) 
1994. 

[8] For a critical discussion, see H. Monien, preprint |cond-mat /0 1 1 1 7*81 

[9] N. Dupuis, Phys. Rev. B, 65 (2002) 245118. 
[10] Z.Y. Weng, C.S. Ting and T.K. Lee, Phys. Rev. B, 43 (1991) 3790. 
[11] F.D.M. Haldane, Phys. Lett, 93A (1983) 464. 

J12] D. Baeriswyl (Editor), H.J. Schulz, The Hubbard Model (Plenum, New York) 1995. 

113] E. Manousakis, Rev. Mod. Phys., 63 (1991) 1. 

J14] L. Bartosch and P. Kopietz, Phys. Rev. B, 62 (2002) R16223. 

[15| S. Moukouri and M. Jarrell, Phys. Rev. Lett. 87, 167010 (2001). 



